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In inflationary cosmology, the particles constituting the Universe are created after 
inflation in the process of reheating due to the interaction with the oscillating 
inflaton fleld. We briefly review the basics of the slow reheating, and the stage 
of fast preheating, when the particles are created explosively in the regime of 
parametric resonance. The non-perturbative, out-of-equilibrium character of the 
parametric resonance changes many features of reheating. For these proceedings, 
we will highlight a few aspects of preheating: the structural dependence of the 
parametric resonance on the inflationary model, including V((f>) = ^^(f)^, ^0*, 
1 — cos ^; "rescattering" of created particles; and phase transitions after inflation. 



1 Reheating after Inflation 

In modern versions of inflationary cosmology there is no pre-inflationary hot 
stage, the Universe initially expands quasi-exponentially in a vacuum-like state 
without temperature. During inflation, all energy is contained in the inflaton 
field (j) which is slowly rolling down to the minimum of its effective potential 
V{4)). When chaotic inflation ends dX (f) ^ Alp, the inflaton field begins to 
oscillate near the minimum of V{(j)) with a very large amplitude 4> ~ j^Mp. 
In this scenario all the particles constituting the Universe are created due to 
the interaction with the oscillating inflaton field. Gradually, the energy of 
inflaton oscillations is transferred into energy of the ultra-relativistic particles. 
Eventually created particles come to a state of thermal equilibrium at some 
temperature Tr, which is called the reheating temperature. The transitional 
stage between inflation and the standard Big Bang is called the period of 
reheating after inflation. 

To describe creation of elementary particles from the inflaton oscillations, 
we shall consider the interaction terms in the fundamental Lagrangian. The 
inflaton field (f) may decay into bosons x fcrmions ip due to the interac- 
tion terms — ^g^^^X^ ^^'^ —h'lljipcj), or into its own Bose quanta 6(j) due to the 
self-interaction X4P'6(lP , or due to the gravitational interaction; A, g and h are 
small coupling constants. In the case of spontaneous symmetry breaking at a 
scale (T, the term —^g^cji^x^ gives rise to the three-legs term —g'^a(f>x'^- We 
will assume for simplicity that the bare masses of the fields x ^-nd ip are very 
small. The elementary theory of reheatiag is based on the perturbative theory 
with respect to the coupling constantsll^. Inflaton oscillations are interpreted 
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as a superposition of a number of (/)-particles, inflatons, at rest. Each infla- 
ton has an energy equal to the frequency of the background osciUations. The 
decay rate of inflatons (or the particle production rate) is given by perturba- 
tion theory. For simplicity, consider the interaction —^g^F{(j))x^ between the 
classical inflaton field (j) f^nd the quantum Bose field x with the eigenfunctions 
Xk{t) e~^^^ with comoving momenta k. The temporal part of the eigenfunction 
obeys the equation 

Xk + 3-Xk+(^+9^F{c^))xk^0 , (1) 



a \ a 



with the vacuum-like initial condition Xk — ^yj^ the far past. The WKB 



solution of Eq. (|l|) is 



jdt 
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where the time-dependent frequency is Lo1.{t) = ^ + g^c})^ plus a small correc- 
tion ^ H^,H; initially f3k{t) = 0. For \f3k\ <C 1, an iterative solution 

t t' 
(3k — \ J dt' exp (— 2i J dt"Lj{t")^ ^ gives the standard result of per- 

— oo — oo 

turbation theory for the particle occupation number nj- = |/3fcp. The integral 
can be evaluated by the stationary phase method. For the three-legs interac- 
tion g'^uipx^) the perturbative result can be interpreted as separate inflatons 
decaying independently of each other into pairs of x-particles. The resulting 
rates of the three-legs processes cf) — > XX ^^id <f) — > plus gravitational decay 
are given by 

- 8^^^ ' ^-V-V - 8^ - i 9 - 8^^2 ■ (■^) 

For the four-legs interaction a pair of inflatons is decaying independently into a 
pair of x-particles, {(fxj) —^ xx)- However, the decay of massive inflatons in this 
case is different. The transfer of energy from inflatons to the created particles 
rapidly decreases with the expansion of the Universe as ^^(q^^x) ^ a~^. 
Therefore, the complete decay of the massive inflaton field in the theory with 
no spontaneous symmetry breaking or with no interactions with fermions is 
impossible. 

Reheating completes at the moment tr when the Hubble parameter H — ^ 
becomes equal to Ttot- Assuming thermodynamic equilibrium sets in quickly 
at the temperature T^, one can equalize the energy density of the universe and 
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thermal energy of created ultrarelativistic particles at the moment tr- This 
gives the reheating temperature 

Tr ~ O.lv/TMp , (4) 

which is ultimately determined by the total decay rate Ttotai- In order to 
get numerical estimates for tr, Ttot and Tr, one should know the frequency of 
the inflaton oscillations and the coupling constants. The coupling constants 
of interaction of the inflaton field with matter cannot be too large, otherwise 
the radiative corrections alter the shape of the inflaton potential (unless SUSY 
eliminates radiative corrections). The necessary condition for the decay of 
inflaton oscillations is that their frequency is greater than the effective masses 
of created particles. Parameters of the inflaton potential are restricted from the 
constraints on amplitude of the cosmological fluctuations. All together it allows 
us to put constraints on the perturbative reheating. The largest total decay 
rate is F < 10~^*'Mp. It takes at least 10^'' oscillations to transfer the energy of 
inflaton oscillations into the created particles, i.e. this reheating is very slow. 
The bound on the slow reheating temperature is Tr < 10^ GeV . This is a very 
small temperature, at which neither the standard mechanisms of baryogenesis 
in the GUTs work, nor cosmologically interesting heavy strings, monopoles and 
textures can be produced. In the slow reheating scenario with the potential 

2,2 

, the post-inflationary stage of inflaton oscillations with a scalar factor 
a{t) ^ t^/^ lasts sufficiently long. It can enhance the gravitational instability 
of density fluctuations, which leads to the formation of the primordial black 
holes for some specific spectra of the initial density fluctuations. 

2 The Stage of Preheating 

In the elementary theory outlined above, we made signiflcant oversimplifica- 
tions, assuming that each inflaton decays independently. However, interacting 
with quantum particles, inflatons act not separately jbut as the coherently oscil- 
lating homogeneous field as follows from Eq. The oscillating effective 
frequency LLi{t) and very large amplitude of the background oscillations (t){t) 
can result in a broad parametric resonance of the "oscillator" amplitude Xk- 
The smallness of alone does not necessarily correspond to small occupation 
number \[3k\, which can grow exponentially. Then the energy is very rapidly 
transferred from the inflaton field to other bose fields. This process occiixs 
far away from thermal equilibrium, and therefore it is called preheatingux^sa. 
Reheating never completes at the first stage of parametric resonance; eventu- 
ally the resonance becomes narrow and inefficient. The inflaton fleld decay is 
completed and created particles are settled in thermal equilibrium during the 
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subsequent, final stages of reheating. However, the elementary theory should 
be applied not to the original inflaton oscillations, but to the products of its 
decay, or to its residual oscillations. The short stage of explosive preheating in 
the broad resonance out-of-equilibrium regime may have long-lasting effects on 
the subsequent evolution of the universe (compare to other cosmological epochs 
where particles temporarily can be out-of-equilibrium). It may lead to specific 
nonthermal phase transitions in the early universe and to the production of 
toQQlQgical defects □, it may make possible novel mechanisms of baryogene- 
sisEj'p, and it may change the final value of the reheating temperature T^. 

The non-perturbative character of the parametric resonance changes many 
features of the particle creation from inflatons. For these proceedings, I will 
highlight only several interesting aspects of preheating, which we have learned 
recently: the structural dependence of the parametric resonance on the theory 
of inflation V {(/)); "rescattering" (mode mixing) effect for created particles; and 
phase transitions after inflation. 

3 Structure of Resonance and Inflaton Potential V{(p). 

The resonance depends on the form of the background oscillations, which are 
defined by V{(f)). We consider several popular models of ¥{(/)). 

3.1 Resonance in Inflation 

For the quadratic inflaton potential we have the background inflaton oscil- 
lations (j){t) ~ 0o(i) • sin(m,/,i), (j)o{t) = ■ The scalar factor is 

a{t) « aot^^^- The amplitude of oscillations is decreasing as the universe ex- 
pands, (t>Q{t) ^ O.lA/nfl"'^/^ Let us consider ^(jj^x^ interaction. The equation 
for the fluctuations (|^) with F{(t>) = 0^ in a new time variable z = rat takes the 
form reminiscent of the Mathieu equation. The dimensionless coupling param- 
eter q — depends on time via <i>^ cx a^^. Without the expansion of the 
universe, that equation is exactly the Mathieu equation. Since ~ 10~^, 

we can have q ~ lO^'^^f;^ ^ 1, so the parametric resonance is broad. The 
regular Mathieu resonance can be described by its stability/instability chart. 
The Mathieu equation admits exponentially unstable solutions Xk c< exp(/ij.z) 
within the set of resonance bands, corresponding to the parametric resonance. 
For given parameter g, there are separated stability and instability zones. For 
g 3> 1 the parametric resonance is broad, with the width of the leading reso- 
nance band ~ OA^/gmM and with the characteristic exponent /i^ ranging from 
zero to its maximum 0.28. In this case the eigenfunctions Xfe(t) have the form 
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Figure 1: Left: characteristic exponent of the Mathieu equation as function of the 
parameter q for the fluctuation with fc = 0. Right: the characteristic exponent fij_ in an 

expanding universe as function of k for j = 5 and 6. fij. changes dramatically with j but 
within the theoretical envelope 



of the WKB solution (|l|) for all t except very short time around moments tj, 
j = 1,2,3,..., where the background oscillating field crosses zero, = 0. 

Non-adiabatic changes of the amplitude Xk (t) occur only in the vicinity of tj , 
which can be described with the parabolic scattering. An essential parame- 
ter here is the phase 9k = Jf'*^ dt^Jk'^ + g'^(j)'^{t) accumulating between two 
successive zeros of (/>(i). If the parameter q = const, the phase 6^ as function 
of k is not varying. Then we have a stable phase correlation/anticorrelation 
between successive scatterings at the parabolic potentials, which produces the 
sequence of stability/instability bands, see Fig. |^. 

However, in an expanding universe the parameter q is time-depended, q = 
^J^a ■ Q oc For the broad resonance case q ^ I this parameter can jump 
over a number of instability bands within a single oscillation, and the concept 
of stability/instability bands is inapplicable here. It is easy to understand that 
in this case the parametric resonance is a stochastic process. Indeed, for large 
initial values of q, the phase rotation between tj and ij+i, S9k — ^ is much 
larger than tt and therefore can be considered as a random nuniLer in the 
interval [0, 27r]. As a result, there are no separate stability/instability bands, 
but each mode can be amplified or deamplified with every half period of the 
infiaton oscillation. An example of a stochastic fik is plotted in Fig. |^. Positive 
and negative occurrences of Hk for k < Tr~^ appear in the proportion 3:1, and 
the net effect is parametric amplification. 
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3.2 Resonance in jX(f>'^ Inflation 

For the theory with the potential V{4>) — j^^ip"^ it is more convenient to express 
the background solutions via the conformal time t = J and conformal field 

ip — acf) variables: (^(t) k (p cn ^r, , a(r) = sJ^-^t, (p is the amplitude 
of oscillations. The oscillations in this theory are not sinusoidal, but given by 
elliptic function. Let us consider a ^(f^y^ interaction of the inflaton. Eq. 
for quantum fluctuations Xk can be simplified in this theory. Using a conformal 
transformation of the mode fimction Xk{t) = a{t)xk{t), from Eq. (|^) we obtain 



A V ' 

where = and g = ^. The equation for fluctuations does not de- 

pend on the expansion of the universe and is completely reduced to the similar 
problem in Minkowski space-time. This is a special feature of the conformal 
theory iA^l^dr \g'^^x^ ■ The mode equation (||) belongs to the class of Lame 
equations Q'DO The combination of parameters q = jX ultimately defines 
the structure of the parametric resonance in this theory. This means that 
the condition of a broad parametric resonance does not require a large initial 
amplitude of the inflaton field, as for the quadratic potential. The strength 
of the resonance depends non-monotonically on this parameter. The stabil- 
ity/instability chart of the Lame equation (H) in the variables ^k^,^ 

constructed inS (see also P.Greene's contribution to this volume). On the 
left panel of Fig. H we show slices of the Lame stability/instability chart as a 
functions /Xfe for several values of parameter 

To see how the general theory of the Lame resonance works, for illustration 
let us consider parametric resonance in the model of the self-interacting two- 
component scalar field (</ii, </i2) with the effective potential V{<^) = ■|(0i+</'2)^- 
For the two component inflaton scalar field, we always can align the classical 
background field with one of the components, say 4>i{t). The equations for 
the quantum fluctuations for two components S4>i and 6(j)2 are different. The 
equation for the mode function of fluctuations in the direction 0i is 

(/-ifc+a- </)ife+ ^ + 3A</)? Uifc = . (6) 

The equation for the mode function of fluctuations in the direction 02 is 

+ hk+{^ + X(j)l](j)2k=0 . (7) 
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Figure 2: Left: The characteristic exponent for the Lame equation ^ as a function 
of (in units of Xip^ ) for ^ = 1.0,1.5,2.0,2.5 and 3.0, labeled a through e respectively. 
Right: fii^ as a function of (in units of m?) in the Sine-Gordon theory. Different curves 
correspond to different values of the background misalignment angle, 9o, from the outer 
curve inwards 9o = 1-67, 1.25, 1.0, 0.83, 0.71, 0.625, 0.56, 0.5, 0.45, 0.42, 0.38, 0.36. 



In some early papers on preheating the iteration series in the small parameter 
A was used. This approach is incorrect. Indeed, using the conformal transfor- 
mation of the mode functions ipik{t) = a{t)(j)ik{t) and (p2k{t) = cL{t)(j)2k{t) , we 
can reduce the equations (|6|), to the general equation with particular 
value of the parameter q = 3 and 1 respectively, so there is no small parameter 
associated with the resonance. In both cases q = 1,3 thcjcje-is only one insta- 
bility band, but the strength of the resonance is different ala^cl. The resonance 
in the "inflaton" direction 0i is weak, maximal value of the characteristic ex- 
ponent of the fluctuations (pik oc e^^'^ is /ii « 0.036; the resonance in the other 
direction (/)2 is much stronger and broader, /i2 ~ 0.147, see Fig. |^. 



3.3 Resonance in Sine- Gordon theory 

The Sine-Gordon potential V{(i>) = m^f^ ^1 — cos is the prototype of the 
natural inflation scenario, as well as the cosmic axions. Let us consider small 
quantum fluctuations Scj) around oscillating homogeneous background field 
(j){t), which can be excited due to the self- interaction. At the moment we 
neglect the expansion of the universe. The equation for the temporal part of 
the eigenfunctions of fluctuations (j)k {t) is 

(fk+ (^k^ +ml cos j^(j)k=0 . (8) 
7 



Since the background oscillations jare periodic, there is parametric amplifica- 
tion of the fluctuations (jik (t) . We □ have derived the general analytic solution 
of this equation for an arbitrary initial amplitude 0o, or for an arbitrary mis- 
alignment angle 9 = The characteristic exponent fi^ as function of k is 
plotted in Fig (^). There is always only a single resonance band. The value 
of fj,k strongly depends on the phase Oq. In the important case of small mis- 
alignment angle Oq ^ 1, the general formula for /i^ is reduced to the simple 
expression /ife = j^yj&o ^ In this case the maximum value — ^ 

— The resonance band in this limit is < < It is instructive 
to return to the equation for fluctuations (|^), and perform there a small 9 
approximation to the cos term. This gives us the Mathieu equation with the 
parameter q = We would predict a set of resonant bands, the first reso- 
nance band being located at -s^- < < which is shifted from the actual 

^ 8 — m-^ — 8 ' 

instability zone. The difference between the actual Sine-Gordon resonance and 
its Mathieu approximation is even more noticeable if we include the expansion 
of the universe. The parametric resonance due to the self-interaction is not 
effective both for the natural infiation and for axions, because the amplitude 
is quickly 6*0 ~ a~'^/^ decreasing with the expansion, but not because of the 
redshift from the resonance as was thought previously. Certainly, if we allow 
coupling to other bosons, say a (;^(/)^x^"iiiteraction, then preheating in natural 
inflation can be very efficient. 

4 "Rescattering" of created particles 

So far we have treated the fluctuations x of S4> as test fields in a given back- 
ground (l){t) and a{t). Due to the exponential instability of fluctuations, one 
expects their backreaction on the background dynamics is gradually accumu- 
lating until it changes the resonance itself. There are two especially important 
effects of backreaction. First, fluctuations may change the frequency of back- 
ground oscillations, which can be taken into account in the Hartree approx- 
imation tM. The second effect is the production of the inflaton fluctuations, 
which occurs due to the interaction of created particles x (or Scj)) with the in- 
fiatons at rest. One can visualize this process as scattering of x particles on the 
oscillating field taking infiatons away from the homogeneous oscillating 
condensate. Fluctuations of Bose fields generated with large occupation num- 
ber can be considered as classical waves with gaussian statistics. Therefore, all 
the field Xi ^4' can be treated as interacting classical waves. Thk-ppens 
an avenue to study preheating through lattice numerical simulations bI13. 
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The theory of "rescattering" can be iUustrated with the two-component 
self-interacting model ¥{(/}) = j{4>i + 4'2)'^ i where the inflaton field is identified 
with the first component 4>i . In this model we have the cross-interaction term 
A(/)f0|. The rescattering of the 02 particles on the classical field (pi leads to 
the production of (pi particles in the process ^2 fc '^i ~* ^2 P gen- 
eration of fluctuations (pi due to this process can be incorporated in the mode 
equation (0) with the scattering term X(Pi{t) J (fik'cp.^ k-k''t'2 k'- conformal 
variables, the fluctuations ipik due to the rescattering are given by the formula 
(cf.li): 

^iAr) = -j7r^ r dr' smcjk{T-T')MT')d'k'^^^_p{T')v^^,{T') , (9) 

[^Tt) OJk Jq ' ' 

where ujk is the effective frequency of Eq. (^). From this it follows that the 
amplitude </3i,fc(T) grows with time as e^^^'^, where ^2 = 0.147. Therefore 
(ipf) evolves in time even faster than ((^2)- Analysis of Eq. (j^) shows that the 
process 6<p^ k 'Pi ^ ^<p2 p ^'P'l k~k" strictly speaking, does not correspond to 
the "particle-like" rescattering, but rather to the interaction of waves, mode 
mixing. 

5 Phase Transitions from Preheating 

The theory of cosmological phase transitions is one of the main topics in mod- 
ern cosmology. Cosmological phase transitions in GUTs, which may happen at 
Tc ~ 10^^ — 10^^ GeV, could give rise to primordial monopoles and other topo- 
logical defects, which could either kill popular cosmological models, or help 
them. The pcia^ibility of the phase transitions from preheating is based on the 
following ideau. The number of particles created during preheating can be very 
large. These particles initially are far away from thermal equilibrium. They 
may change the shape of the effective potential, which may lead to specific 
nonthermal phase transitions soon after inflation. Nonthermal phase transi- 
tions induced by preheating are similar to the usual high-temperature phase 
transitions, but they may occur even in the models where the high tempera- 
ture effects cannot lead to symmetry restoration. This is because fluctuation 
amplitudes (<p'^) from resonant amplification can be very large. 

As an example, let us consider a two-component scalar field ((pi, (P2) with 
the effective potential 

n0) = ^('^?+02-vY. (10) 

If we identify the (pi component with the "inflaton" direction (as in section 3.2), 
the effective mass of the background field is m'^ff = — v^ -I- 3X(pl + 3A((50^) -|- 
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\{5(j)\). We shall consider the growth of fluctuations in this model. At the be- 
ginning one can neglect the tachyonic term. Then we can use the results of Sec. 
3.2, {5(j)\) oc e^^^^ , (<5<^2) ^ e^^^^ ^ ^ Mi- However, the mode mixing due to 
the cross-interaction leads to the faster growth of fluctuations oc e'*'^^'^. 

If > 10~'*i^^, the generation of fluctuations will be completed before the 
bifurcation in the evolution of the the background field occurs at 4>i ~ v. By 
that moment, we will have the amplitude of fluctuations comparable to that 
of the background field, ~ {S(t>\) ^ 4>i- In this model we have symmetry 

restoration alongside with the residua oscillation of the background compo- 
nent (j)i{t), and what is most important, the formation of strings. Notice that 
the string formation mechanism is somewhat different from the Kibble mech- 
anism. Lattice simulations of the self-consistent non-linear disuamics clearly 
demonstrate the formation of topological defects in the model tHi. 
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